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Abstract 

Many microorganisms swim in a highly heterogeneous environment with obstacles such as fibers or poly¬ 
mers. To better understand how this environment affects microorganism swimming, we study propulsion of 
a cylinder or filament in a fluid with a sparse, stationary network of obstructions modeled by the Brinkman 
equation. The mathematical analysis of swimming speeds is investigated by studying an infinite-length 
cylinder propagating lateral or spiral displacement waves. For fixed bending kinematics, we find that swim¬ 
ming speeds are enhanced due to the added resistance from the fibers. In addition, we examine the work 
and the torque exerted on the cylinder in relation to the resistance. The solutions for the torque, swimming 
speed, and work of an infinite-length cylinder in a Stokesian fluid are recovered as the resistance is reduced 
to zero. Finally, we compare the asymptotic solutions with the numerical results obtained from the Method 
of Regularized Brinkmanlets. The swimming speed of a finite-length filament decreases as its length de¬ 
creases and planar bending induces an angular velocity that increases linearly with added resistance. The 
comparisons between the asymptotic analysis and computation give insight on the effect of the length of the 
filament, the permeability, and the thickness of the cylinder in terms of the overall performance of planar 
and helical swimmers. 
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I. INTRODUCTION 


The self-propulsion of microorganisms that utilize flagellar propulsion has been the topic of 
a vast number of analytical, experimental, and computational studies for many years (reviewed 
in p). Many species of spermatozoa and bacteria are able to swim by propagating lateral or 
spiral waves along their cylindrical flagella [2]-4]. Similarly, larger organisms such as C. elegans 
(nematodes) are also able to make forward progression through soil via undulatory locomotion [[5]. 
The native environment in which these organisms live varies greatly. For example, spermatozoa 
encounter different fluid environments in the female reproductive tract that include swimming 
through or around mucus, cells, hormones, and other large proteins [J6HB3I- Similarly, bacteria are 
able to swim in the mucus layer that coats the stomach and move in biofilms with extracellular 
polymeric substances Il2l l9l [TOl. 

One may wonder, how does the swimming speed or mode of swimming change in these differ¬ 
ent environments? Early experiments showed that Leptospira, a slender helical bacterium, is able 
to swim faster in methylcellulose (MC), a gel with chains of long polymers ifTTl . Another study 
showed that swimming speeds of seven different types of bacteria were enhanced in higher viscos¬ 
ity solutions of MC and PVP (polyvinylpyrollidone) [fl2l ; beyond a certain viscosity or polymer 
concentration, this enhancement was no longer observed. Experiments of sperm in MC and PA 
(polyacrylamide) gels showed that swimming speeds, beat frequency, and amplitude of undulation 
vary as the viscosity and concentration of the gels are varied [|T3l 141. C. elegans have also been 
observed to swim faster in polymer networks 0. 

Since the length scale of these swimmers is small, they live in a viscosity dominated environ¬ 
ment where inertia can be neglected. Many studies have focused on analyzing idealized swimmers 
in viscous fluids at zero Reynolds number. Seminal work by GI Taylor examined swimming speeds 
of an infinite sheet in two-dimensions (2D) and an infinite cylinder with circular cross section of 
small radius in three-dimensions (3D), propagating lateral displacement waves 03] [16]. In these 
studies, it was shown that the second order swimming speed scales quadratically with amplitude 
and linearly with frequency for small amplitude bending. This analysis has been extended for 
several different cases including swimming speeds for cylinders with non-circular cross sections 
ifTTIl . as well as improvements to the perturbation series lfl8l . 

Since the fluid that these swimmers are moving through contains different amounts of proteins 
or other structures, more complex fluid models have been proposed and analyzed. For the case of a 
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swimming sheet, studies have looked at the asymptotic swimming speeds in a gel represented as a 
two-phase fluid (elastic polymer network and viscous fluid) where enhancement in propulsion was 
observed for stiff and compressible networks IIT91 . In contrast, a two-fluid model (with intermixed 
fluids) exhibited a decreased swimming speed relative to the case of a fluid with a single viscosity 
in both asymptotics and numerical simulations Il20l . In another model, Magariyama el al. Eli 
looked at a fluid governed by two viscosities using a modified resistive force theory and found 
that there is an enhancement in propulsion efficiency when the viscosity of the polymer solution 
increases and the other fluid viscosity is held constant. Swimming in a shear thinning fluid has 
also been studied; locomotion of finite-length swimmers is enhanced (2D numerical simulations) 
ll22l and infinite undulating sheets exhibit a decrease in swimming speed relative to the Stokes 
case ll23l . 

Since the proteins or polymer chains in gels may cause a fluid to exhibit a nonlinear strain 
response (frequency dependent), viscoelastic fluid models have also been considered. Through 
asymptotic analysis, it has been shown that the swimming speed of infinite sheets and cylinders in 
viscoelastic fluids decreases relative to the speed in a purely viscous fluid |[24l - [26l . Simulations of 
finite-length swimmers in a viscoelastic fluid at zero Reynolds number governed by the Oldroyd 
B equation revealed that enhancement in swimming speeds can be observed when asymmetrical 
beatforms and swimmer elasticity work together [[271 [28 ]. Specifically, increases in swimming 
speeds were observed in a viscoelastic fluid when the beat frequency of the swimmer is on the 
same time scale as the polymer relaxation time [|27l ; when the polymer relaxation time is fast, 
other models may be more appropriate to understand swimming speeds. 

Another approach is to think of the fluid with an embedded polymer network as a porous 
medium. Darcy’s law has been used to describe the fluid flow in porous media, where average 
velocity is proportional to the gradient in pressure. This law is not able to capture contributions of 
the viscous stress tensor and it is only valid on the macroscopic scale where the domain is large and 
boundary effects can be neglected |[29l 30]. To overcome these disadvantages, the incompressible 
Brinkman flow equation has an additional diffusion term 11291 . 

Vp = /iAu — —u, V ■ u = 0, (1) 

7 

where p is the pressure, u is the velocity of the fluid, 7 is the permeability of the porous medium, 
and the effective viscosity is ji. This equation represents the effective flow through a network of 
stationary obstacles with small volume fraction [ 29[ [3H - I331 . The resistance due to the obstacles 
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is characterized by /if 7. Note that the incompressible Stokes equations are recovered in the limit 
as 7 —>■ 00, and when 7 —> 0, Eq. (|TJ) will behave like Darcy’s law. Another characteristic of a 
Brinkman fluid is the Brinkman screening length, ^ 7 , which marks the approximate length over 
which a disturbance to the velocity would decay. For comparison, in 3D, the flow due to a point 
force in Stokes flow decays as 1/r whereas the flow due to a point force in a Brinkman flow decays 
like 7/r 3 D 31 I 351 . 

In the case of a two-phase fluid composed of a polymer network and solvent, if the polymer is 
stationary, we obtain the Brinkman equation. In this limiting case of a two-phase fluid, an infinite- 
length sheet exhibits an enhancement in swimming speed [|T9ll . Previously, Leshansky ll35l derived 
the asymptotic swimming speed for an infinite sheet propagating waves of lateral bending in a 
fluid governed by the Brinkman equation. They observed that swimming speeds scaled similarly 
to those of Stokes, scaling quadratically with amplitude. In addition to the Stokesian swimming 
speed, there is an extra factor that depends on the permeability and is monotonically increasing for 
decreasing permeability (increasing the resistance in the fluid). 

In this paper, we focus on calculating the asymptotic swimming speed for a waving cylindri¬ 
cal tail that exhibits lateral displacement waves in a Brinkman fluid. A second order asymptotic 
swimming speed is derived for planar bending and we find that swimming speeds are enhanced, 
similar to the 2D case for an infinite sheet. Swimming speeds are also calculated for cylindri¬ 
cal tails with spiral displacement waves, showing that fluid resistance enhances swimming speed. 
These results shed insight on how added fluid resistance changes propulsion of cylindrical tails 
when the kinematics are prescribed. In addition, as the resistance approaches zero, we recover the 
swimming speed, work, and torque for an infinite-length cylinder in a fluid governed by the Stokes 
equation. Through our analysis, we also find the range of enhancement in swimming speeds for 
the infinite cylinder in a Brinkman fluid and the relation to permeability, cylinder thickness, and 
wavenumber. To validate our asymptotic results, we apply the method of Regularized Brinkman- 
lets to study finite-length filaments. Through validation, we find that the theoretical swimming 
speed of filaments with planar bending waves matches up well with the simulation data and that 
the asymptotics overestimate swimming speeds for shorter length cylindrical swimmers. In the 
helical bending wave case, we calculate the external torque exerted on the filament by the sur¬ 
rounding fluid. We observe that the numerical and the asymptotic findings may not consistently 
agree with one another; the asymptotics overestimate the torque of finite-length helical swimmers. 
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II. SWIMMING SPEEDS FOR A CYLINDER WITH PLANAR BENDING 



FIG. 1. Current configuration (deformed state) of the cylinder propagating bending waves is shown with the 
solid circle. O' is the origin for the current configuration and O is the center of the original (non-deformed) 
state that is shown with the dashed circle. 


A. Cylinder with lateral displacement waves 

Similar to previous work [TUGS, we consider a cylinder of constant cross section, bending 
with small amplitude in the x direction, immersed in a fluid. The cylinder is bending in the x — y 
direction with 

x = bsm(k(z + Ut)), y = 0, 

where b is the amplitude, U is the velocity of the propagating wave, and k is the wavenumber, 
defined as k = 2tt/A where A is the wavelength. With this, the velocity components of the 
cylinder have the form u x = bkU cos (k(z + Ut )) and u y = 0. To simplify, we let s = k(z + Ut) 
and convert the above equations into cylindrical coordinates to obtain the boundary conditions on 
the surface of the cylinder, 

u r = bkU cos 0 cos s, ug = — bkU sin 6 cos s, u z = 0. (2) 

From this point, we will regard the velocity components in cylindrical coordinates as u r = u, u g = 

v and u z = w. 

The time-dependent position of the cylinder at any given point on the surface is given as 

r 2 = a 2 + b 2 sin 2 s + 2 ab sin s cos 9', (3) 
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as shown in Fig. [l]for 9' = 0 + £. As detailed in Appendix |IXA[ we can show that 


r = a + b sin s cos 9 


( 4 ) 


in the first order of b/a. 


B. Fluid Model 


The 3D Brinkman equation in cylindrical coordinates is: 
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(5) 

( 6 ) 

(V) 


where u, v and w are the velocity components in the direction of r, 9, and z, respectively. The 
continuity equation for the incompressible flow is given by 
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( 8 ) 


Taking the divergence of Eq. ([5]) and using Eq. ([ 8 ]) to simplify, we find that the pressure satisfies 
V 2 p = 0. Let Q — hr and recall s = k(z + Ut). The general solution for the pressure is thus 


p = iiknA mn K m (n () cos (m9) cos (ns), (9) 

where K m is the m th order modified Bessel function of the second kind and A mn is a constant 
which is evaluated using the boundary conditions Il36l . Based on the pressure in Eq. (|9]), we 
assume the velocity components can be described as 


u = u mn cos rn9 cos ns, v = v mn sin m9 cos ns, and w = w mn cos m9 sin ns. (10) 


Note that u mn , v rnn and w rnn are functions with respect to ( only. Substituting u, v, w, and p from 
Eqs. into Eqs. ([5])-([6]) and using the relations s = k(z + Ut) and ( = kr, we obtain the 

following system of equations: 
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with a 2 = I/7 (where 7 is the Darcy permeability). The parameter a is known as the hydrody¬ 
namic resistance of the porous medium and has units of inverse length. In addition, a is propor¬ 
tional to the ratio of the diameter of the fiber over the spacing within the network. This ratio is 
usually characterized as the mesh spacing Il34ll . 

The homogeneous solutions for Eqs. ([TT|)-([T2|) include the modified Bessel function of the first 
kind, which will diverge as ( —> 00. Thus, we eliminate this solution to maintain finite values for 
the velocities. The particular solutions are 


A. 


A 


i^mn “ 1 “ Vmn)p ^9 ^ and n mn )p „ (13) 


P 


P 2 


where (3 = a/k is the scaled resistance. It is a nondimensional constant that characterizes the 
relationship between the resistance or average mesh size and the wavelength of the swimmer. 
After simplifying, the general solutions to Eqs. ([TTT)-([T2|) are 


Hmm T Ann 


'U'mn Vmn 


B mn K m+1 ( X C) + -j^n 2 K m+l (n(), 
C mn K m _ i(xC) + 


(14) 

(15) 


for x = s/n 2 + /3 2 . The constants B mn and C mn are determined by the boundary conditions of the 
cylindrical tail. The radial and tangential velocity components are found to satisfy the following 
equations: 

2 U mn = B mn K m+ i(xC) + CmnKm- l(xC) + [Km+lijlC) + K m _i(n()\ , (16) 

P 

2 a 

2v mn = B mn K m+ i(xC) - C mn K m _ i(xC) + noT m nK m {n( ). (17) 

P C 


The axial component of the velocity is determined using the continuity condition given in Eq. ([8]) 
and is given by 


du mn 1 

Ti ll ' n in g T ~ ( M inn T FriV m/n J 

A, p n K m (n() + ^(B mn + C mn )xK m (x(). 


(18) 


Since our goal is to determine the swimming speed of the cylinder, we will have to determine the 
first and second order solutions, using the condition that the disturbance caused by the cylinder 
body should vanish at infinity [fl6l . 
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C. First order solution 


As detailed in Appendix |IXB[ the velocity components are expanded about ( = 
the first order, when rn — 1 and n = 1 , the boundary conditions are u , = bkU, V\ 
wi = 0. Plugging into Eqs. (fl4])-(fT5]) and ( |T8| ), we obtain: 


Ci = ka. To 
= —bkU, and 


u 1 + v 1 = BK 2 (xCi) + ^-MCi) = 0, 

Ul - Vl = CAo(xCi) + ^Ao(Ci) = 2 bkU, 
w i — T^-^i(Ci) + + C)xBi(xCi) — 0, 

for x = y/l + p. From Eqs. ([I9j)-([2T]), the constants are 

A 1 2 bkU 

^ __ 0(Ci)iTo(xCi)’ 

1 2 bkU Ci) 

0(Ci) Ao(xCi) ^CxCi)’ 
r ZbkU r 1 iT 0 (Ci) 1 

■ Ao(xCi) L %(Ci)Ao(xCi)J’ 

where 


(19) 

( 20 ) 
( 21 ) 

( 22 ) 

(23) 

(24) 


2 AC (Ci) Aq(Ci) K 2 ((i) 

xKiixCi) KoixCi) K 2 { X CiY 


(25) 


To determine the velocity of the cylinder, we have that Eqs. (|70|)-(|72]) in Appendix |IX B| will vanish 
at infinity llT6ll . Thus, there is no contribution to the swimming speed of the cylinder in the first 
order expansion. 


D. Second order solution 


The second order expansions and boundary conditions are detailed in Appendix |IX B[ Using 
the same argument for the velocity of the filament at infinity, we arrive at 

t/oo ~^bk(^w^) ^ = ^ 1= i ca , 

where w\ is the first derivative of the axial velocity component given in Eq. ( fl8| ) with respect to C 
(for m — 1, n — 1). Using the first order solution, and evaluating w[ at the boundary, C = Ci = 
we have 


2 bkU r ATq(Ci) , ATi(Ci) ' 
tt(i) L^o(xCi) X Ki(x Cl). ■ 


(26) 
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The swimming speed Uoo up to second order expansion is thus 


= 2 bVu m 


Ao(Ci 


X- 


A\(Ci 


(27) 


L^o(xCi) -MxCO 

The asymptotic velocity for an infinite-length cylinder that is propagating planar bending waves in 


a Brinkman fluid is given above in Eq. ( [27] ) and depends on the scaled resistance a/k through y. 
In the limiting case when Ci = ka <C 1, the limit forms of the Bessel functions are [[371 : 
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Sl 
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*o(Ci) = - lo s Ci - 7e + log 2 + e>(Ci 2 ), 
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where y e is the Euler-Mascheroni constant. 

Thus, for Ci 


xKl 


+ o(xCt), 


log Cl + log 2 - 7 e - logX + 0(x Cl), 


0(Ci) = 2-x 2 


Ao(Ci) 


K 0 (Ci) ~ 1 °SX' 


To second order, the nondimensional swimming speed, U^/U, in the case of a cylinder propagat¬ 
ing lateral bending waves is given as 


Up o _ 1,2.2 [ (l~X 2 )-MCi) +X 2 logX 

U 2 [(1 -x 2 )K 0 {Ci)-(2-X 2 )logx. 


(28) 


for ka <ti 1. We note that this swimming speed scales quadratically with the amplitude of bending 
b and depends on the resistance a through the parameter y. The swimming speeds are shown in 
Fig. [2]for several permeability values 7. For comparison, we also plot the swimming speed of the 
same infinite-length cylinder propagating planar bending in a fluid governed by the incompressible 
Stokes equation, as derived by Taylor [[T6l . We observe in Fig. [2] that as a —>■ 0 (or 7 —* 00), we 
approach the Stokes swimming speed. In the next section, we will study this case further. 


E. Comparison of Swimming Speeds 

The Brinkman equation is characterized by the Darcy permeability 7. In the case of 7 — > 00 (or 
resistance a —> 0), we recover the Stokes equation. To understand what happens to the swimming 
speed of the infinite-length cylinder (with ka <C 1) as a —> 0, we will work with Eq. ([28]) to obtain 


9 















(a) (b) 


FIG. 2. (a) The nondimensional swimming speed of a cylinder with planar undulations, calculated from 
Eq. ( [28] ), is shown for several permeability values 7 for fixed wavelength A = 24 and a = 0.05. The Stokes 
case is also plotted for comparison, (b) The ratio between the swimming speed in Eq. ( [28] ) with the speed 
of a swimming sheet derived from Il35ll for different scaled resistance a/k and with k = 2tt/24. 


the following expression, 


^ = - b 2 k 2 
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We note the following limits as a —* 0: 


(29) 
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lim log 1 + — = 0. 
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(30) 


Thus, the second order asymptotic velocity of a cylinder with £1 = ka 1 in a Brinkman fluid 
becomes 

r 


U °° = 

U 2 b k 


Ao(Ci 


tfo(Ct) + j 


This is precisely the asymptotic velocity of the same cylinder immersed in a fluid governed by the 
Stokes equations as derived by Taylor lfT6l . 
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Next, we study the swimming speed of the infinite-length 3D cylinder in comparison to the 2D 
sheet, where both are propagating planar bending waves. The propulsion of an undulating planar 
sheet was studied by Leshansky 0511 and the swimming speed U Les was found to be: 


^ = irtVi + <* 2 M 


(31) 


for a 2 = I/ 7 . The ratio of Eq. (29) and pT) is 
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1 + ¥ 


We plot U r atio versus the scaled resistance a/kin Fig. [2jb). We observe that the ratio decreases as 
a increases. This implies that the 3D infinite-length cylinder swims slower than the 2D sheet in a 
fluid with the same Darcy permeability. When a/k —>• 0, we see that the ratio approaches 

_ A'q(Ci) - 1/2 
A'o(Ci) + 1/2’ 

for a fixed p. This is the ratio of the swimming speeds of the infinite-length 3D cylinder and 2D 


sheet in a fluid governed by the Stokes equation. 


III. ENERGY TO MAINTAIN PLANAR BENDING 


The force on the surface is calculated as F = cr ■ n where cr is the stress tensor and n is the 

normal vector. The velocity components of u at the boundary r = a are given in Eq. ([2]). The stress 

du f dv v 1 du\ . 

tensor components are given by a rr = —p + 2// —— and a r o = // I —- 1 —— . Evaluating 

or \or r r 06 J 

the stress using the first order velocity solution at the boundary gives 

((T rr )= pk [—cos 9 cos s, 

(cr r -e) c=Cl = pk ( sin 9 cos s, 

V a( > / C=Ci 

where r = (/k. Since we consider a fluid with a low volume fraction of stationary and randomly 
oriented fibers, the total stress applied to the filament is assumed to be entirely due to the fluid and 
not influenced by the fibers. This assumption is valid since the distance between the fibers is large 


compared to the radius of the filament. There is further discussion of this in Section VI 
The rate of work done to maintain planar bending is calculated as follows: 


dW = —F • u = pbk 2 U AK^) cos 2 9 + sin 2 9 


cos 2 s. 


(33) 
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Using Eq. (|T7|). the derivative of v\ is: 


a £ = l(-B + C) X K l(x C)-I 


f 2 K 2 (Q + BK 2 ( X C) 


(34) 


where —,B,C are from Eqs. (|22|)-(|24|). The mean value of the rate of work to maintain the 

Q2 


F- 
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f 2bkU 
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For a cylinder immersed in a Brinkman fluid, the mean value of the total rate of work per unit 
length (A) along the surface of the cylinder (r = a) is then calculated as 

finb 2 k 2 U 2 


W = 


k 2 


A ' o(Cl) + 2 l 55 — 1 j l0g l 1 + F 


cr 


(35) 


Ao(Ci 


when is small, and \ — \/l + ff 2 . When the permeability 


where </>(Ci) = 2 - x 2 - „ , . . 

-^o(xCi) 

approaches infinity, the Brinkman fluid behaves like Stokes flow. Thus, when 7 —» 00 (or a —> 0) 


and using Eq. (f30b, we have 


— _ fiTcb 2 k 2 U 2 


Ao(Ci) + l/2' 

This is exactly the same energy contribution to maintain the flagellum in motion in Stokes flow 
m. The nondimensional rate of work is shown in Fig. [3]for several different permeabilities 7 and 
we observe that as 7 gets large, it approaches the work done in a Stokesian fluid. In this analysis, as 
the permeability decreases, we observe that there are small changes in the swimming speed (shown 
in Fig. [IJa)), but the work done increases greatly (shown in Fig. [3]). The mathematical analysis 


for this observation is detailed in Appendix IX D The physical meaning of this phenomenon can 
be explained as follows. For a small permeability, there is a large added resistance present in the 
fluid, preventing the swimmer from propelling itself forward. Therefore, it requires more work to 
move with the same prescribed kinematics. We note that the rate of work of the swimming sheet 
has been previously calculated and is also an increasing function of resistance a 


IV. CYLINDER WITH SPIRAL BENDING 

Next, we consider an infinite-length cylinder propagating spiral waves, motivated by exper¬ 
iments where sperm flagella are able to exhibit helical bending 0. Thus, it is compelling to 
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b 


FIG. 3. The nondimensional rate of mean work done to maintain planar bending along the infinite-length 
cylinder, calculated from Eq. ( [33] ) for several permeabilities 7 where A = 24 and a = 0.05. The Stokes 
case is also plotted for comparison. 


consider the rotational movements of a cylinder propagating spiral bending waves (helical bend¬ 
ing waves with constant radius). One can verify from Fig. [4] that similar to the planar case, to the 
first order of 6/a, 


r = a + 6cos(0 — s), 


(36) 


where details can be found in Appendix IXC and s = k(z + Ut). Eq. ( |36| ) corresponds to a 
cylinder that will achieve the form of a right-handed helix about its axis with angular velocity kU 
in the direction of increasing 6. The formulation for the cylinder is 


x = 6coss, y — bsins, z = bs, 


and the velocity components become 


u x = —bkU sins, u y = bkU cos s, u z = 0. 


Converting the above equations to cylindrical coordinates, we have 


u = bkU sin(6 ) — s), v — bkU co s(9 — s), w — 0. 


(37) 


The motion of the helix includes the contributions of two orthogonal planar motions that are 
perpendicular to the 2 -axis, namely the .x’ 2 -planc and 1 / 2 -plane. The analysis for each plane pro¬ 


ceeds in a similar fashion to that of the planar case, satisfying the boundary conditions in Eq. (37). 
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FIG. 4. Model geometry for a cylinder propagating spiral bending waves. The circle is the cross section of 
the deformed cylinder (current configuration) which is centered at O'. The undeformed cylinder is centered 
at O. 


As previous analysis has shown, the second-order solution can only be determined through first- 
order expansions (see Ifl6l fl71). The second-order velocity components at the boundary are 


U 22 = -^K)?=Ci; v 22 =-\bk(v[) (=(l , W 22 = ^^K)c=d- 


(38) 


Let V s be the propulsion velocity of the helix in the opposite direction of the propagating spiral 
bending waves. With this, similar to Ifl6l . we have 

v * = \ bk K) C =Ci 


where w\ is the same as in Eq. (26). By a simple calculation, we observe 


V s , 2 , 2 Ao(Ci) - X 2 A 0 (xCi) 

U ' (2 -x 2 )K 0 {xCi)-K 0 {CiY 


(39) 


Similar to the results obtained in the planar case, when a —» 0, we recover the speed V s in the 
incompressible Stokes equations, 


K , 2 , 2 Ao(Ci) - 1/2 

U K 0 (( 0 + 1/2' 


Thus, the swimming speed of a spiral bending wave is double that of a planar bending wave with 
the same kinematics. We note that modified resistive force theory calculations have also been used 
to determine expressions for the swimming speed of a spiral bending wave in a Brinkman fluid 

m. 
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In addition to determining the asymptotic swimming speed from spiral bending, we can find 
the expression for the torque exerted on the cylinder by the surrounding fluid. Since the fluid in 
this case flows in a circular motion, the radial and axial velocity components are zero and only 
tangential velocity plays a role in this calculation. That is, 


u — 0, 



Vtk 

T’ 


w — 0, 


where ft is the angular velocity of the helix. With this, we simplify the expression for mean torque 
per unit length applied on the filament by the fluid, = 2irr 2 a r o, to 


Too = —47T/if2. 


(40) 


To solve for Q, we use the boundary condition for n 2 2 in Eq. (38) to obtain 


Vtk 

CT 


-2 6fc ( v 'i)c=Ci- 


(41) 


Substituting Eq. ( [41] ) into Eq. (40) and using Eq. (34) for v\ and Eqs. (22)-(24) to simplify, we 
have 


T = 

1 fYl 


Annb 2 kU 


MCi) + \ - 1 ) i°g 


cr 


of 
+ k 2 


(42) 


In the limit as a —» 0, the torque exerted on the cylinder reduces to 

AirijLb 2 kU 

°° = Ao(G) + l/2’ 

which is the same torque calculated for the Stokes regime by Drummond [|38ll . Note that this 
derivation differs from the work of Taylor llT6l (where w\ was used instead of v \). 


V. RANGE OF PARAMETERS THAT LEAD TO SWIMMING SPEED ENHANCEMENT 


To identify the range of parameter values that lead to enhancement in swimming speeds of the 
infinite-length cylinder with planar waves, we rearrange Eq. (|29]) as follows: 


U 2 tfo(Ci) +1 


1 + 


1 + F -^(CO^iog (f + F 


^o(Ci) - \ 


+ log ( l + i 


(43) 
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Note that Eq. (43) illustrates the velocity behavior in the spiral bending wave case when the 
constant 1/2 is removed. The swimming speed is increasing when the following inequalities hold: 


(i = ka < — exp 

g7e 


la 2 f a 2 

T r i ^ \ 2F og V + F 

Ao(Cl) > a 2 / a 2 

¥- log { 1 + ¥ 


(44) 


la 2 ( a 2 

2 ~k? lQg \ + P 
c? ( a 2 




(45) 


Therefore, for any fixed permeability, the swimming speed of the cylinder is enhanced if the 
thickness and wavelength of the cylinder satisfies the inequality in (|45]). In Fig. [5] we plot the 



FIG. 5. The plot of the function on the right hand side of Eq. (451. 


right hand side of Eq. (45), h(a/k), to show that it is, in fact, decreasing in a manner that is 
dependent on the scaled resistance. This means that if the permeability is reduced, then ka must 
also be reduced to observe swimming enhancement in a Brinkman fluid. Hence, the cylinder 
radius and/or wavelength must decrease in order to observe an increase in swimming speed. This 
finding makes sense since the mesh size decreases as the permeability decreases, thus there is less 
room for the swimmer to move. We note that in addition to an enhancement in swimming speed, 
an increase in torque and rate of work will also be observed when Eq. d45l) is satisfied. 
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VI. RANGE OF PERMEABILITY AND SWIMMING ENHANCEMENT 


Our assumption is that the effective fluid environment can be modeled as a viscous fluid moving 
through a porous, static network of fibers via the Brinkman equation. For small volume fractions 
of fibers, this assumption is thought to be a valid one OH . Further, for randomly oriented fibers, 
Spielman and Goren P3l have derived a relationship between the volume fraction 0, the perme¬ 
ability 7 , and the radius of the fiber a/, as 


Ui r 

-f = 40 

7 


15 a/ Ki(a,f/ ^/y) 
.3 7 GyfyK 0 (af/y/i) 


(46) 


Since the Brinkman model assumes that the fiber network is static, we must have that the distance 
between the fibers (or the interfiber spacing) is large enough for the swimmer to move through 
with little or no interaction with the fibers. To estimate the ratio of interfiber spacing and the fiber 
radius, we use the following equation 051 : 


D 

a f 


2 


/ 1 /3tF 



(47) 


where D is interfiber spacing. In the case where this ratio is large, there are little or no interactions 
between a stationary network of fibers and the swimmers. Thus, it is assumed that the fibers do 
not impart any additional stress onto the filament. 

In Table |T[ we report a few parameter ranges in which we see enhancement of swimming speed. 
In particular, we report ranges of the cylinder radius a, with a fixed wavelength of A = 25 ^m. To 
find these ranges, we use fiber volume fractions and radii from the literature 091 . together with 
our own computed values of permeability from Eq. (46) and average separation from Eq. (|47]). 


Media 

0 

af (nm) 

D (nm) 

7 (f^m 2 ) 

Eq. ( 

45 

), A = 25qm 

Collagen gels, ll39l 

Cervical mucus, 1391 

0.00074 

0.015 

75 

15 

8314 

346 

8.6 

0.0085 

a < 1.337 (nm) 

a < 0.102 (nm) 


TABLE I. The calculated permeability 7 using the given volume fractions 0 and fiber radii a/. The range 
of cylinder thickness a where an enhancement in swimming speed is observed is reported. 


The radii of the principal piece of human, bull, and ram sperm are 0.5 //m. 0.29 //m and 
0.15 //m. respectively Il40l442l| . We note that that flagellar radius decreases along the length of 
the flagellum from the principal piece (closer to cell body) to the endpiece. Thus, swimmers will 
experience enhancement when placed in a collagen gel. However, there will be no enhancement 


17 

















for the three swimmers when they are put in cervical mucus at a volume fraction of 0 = 0.015. 
Further, it is well known that the composition of the cervical and vaginal fluid varies greatly 
through the menstrual or oestrous cycle E|431, and this experimental value of 0 = 0.015 is taken 
at one time point in the cycle lf39ll . For instance, around the time of ovulation, the interfiber spac¬ 
ing D may reach up to 25 /im 0. Using this interfiber spacing and a given fiber radius a/ = 15 


nm, we can further estimate the volume fraction 0 = 0.00033 from Eq. (J47J) and the permeabil¬ 
ity 7 = 0.868 (/im 2 ) from Eq. ( |46| ). Then, the cylinder radii for which enhancement is seen is 
a < 0.765 /im when the wavelength A is 25 /im. At this volume fraction, all three spermatozoa 
species will experience an enhancement in swimming speed in cervical fluid. 


VII. NUMERICAL STUDIES 


A. Background 


In this section, we verify our asymptotic solutions and explore aspects of finite-length swim¬ 
mers using the Method of Regularized Brinkmanlets (MRB) Ii44l . This method is an extensiton 
of the Method of Regularized Stokeslets developed by Cortez [45 , 461 for use with the Stokes 
equations. The general idea is to compute regularized fundamental solutions by replacing singular 
point forces with a smooth approximation. With this, the resulting equations can be solved exactly 
to obtain non-singular fundamental solutions. The smooth approximations to a delta distribution, 
often called ‘blob’ functions, are characterized by a small parameter e that controls their width. 
The singular solutions can be recovered by letting e —>• 0. 

The Brin km anlet is the fundamental solution to the singularly forced Brinkman equation 




Vp = pAu — —u + f 0 <5(x — x 0 ), 


7 


(48) 


where x is any point in the fluid, x 0 is the point where the force is applied, and ()(■) is the delta 
distribution. The pressure and velocity are in the form Il44ll : 


p = f 0 ■ VG(x- xo), (49) 

fiu — f 0 • VVB(x — x 0 ) — f 0 AB(x — x 0 ), (50) 

where G(r) is the Green’s function and B{r) is related to G(r) by the non-homogeneous 
Helmholtz differential equation (A — a 2 )B(r ) = G(r) for r = ||x — x 0 || and a 2 = l/y. 
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(51) 


The solutions of G(r) and B(r) are well known 

G(r) = 1 


47TT ’ 


B(r) = 


1 — e 


4:7ia 2 r 1 


and, thus, the Brinkmanlet velocity in (50) becomes 

fJt u(x) = fo-ffi(r) + (f 0 • (x - x 0 ))(x - x 0 )H 2 (r), 


(52) 


where H\ (r) and // 2 (r) are functions of G(r), B(r), and their derivatives. To regularize the 
fundamental solution, the expression for B(r) is rewritten as 


B e (r) = 




4na 2 R ’ 

where R 2 = r 2 + e 2 so that the singularity is removed. From [f44H . the regularized Brinkmanlet 
velocity is 


where 


A*u(x) = f 0 #i(r) + (f 0 • (x - x 0 ))(x - x 0 )fT|(r), 


g e ~aR 

^ Ana 2 R 5 AixR ! 3 l a 3 R 3 aR ~*” 


H\{(r) = 


D —aR 


+ 


+ —^ + 1 ) + £ 2 Hf, (r 


(53) 

(54) 

(55) 


4vr a 2 R 3 4nR \a 2 R 2 aR 

We note that in the case where the fluid flow is generated due to N point forces, the linearity of 
the Brinkman equation allows the resulting flow to be written as 


N 


/m(x) = M £ (x - x k )f k , 


(56) 


k=l 


where k = 1, • ■ ■ , N and M z [%.) = Hi I + x/,x/, II 2 for x/,. = x — x/,. and identity matrix I. 


Note that x = (x, y, z) is a point in the fluid and force f fc is located at x fc . Eq. (56) is compactly 
written from Eq. ( [53] ) and determines the velocity field on the fluid domain at any given point x. 
Explicitly, f fc = ifl, fjk f^) where the force components are the forces in the x, y and z directions, 
respectively. 


B. Test Cases 

For all test cases, the number of discretization points N depends on the length L of the swimmer 
and are evenly spaced with As = 0.01. Unless otherwise stated, the regularization parameter £ is 
0 . 01 . 
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1. Planar Bending 


We first compare the numerical data obtained from the MRB with the asymptotic swimming 
speed for the case of planar bending. Consider an undulating filament parametrized by the follow¬ 
ing space curve equation as 


x(s,t) — s, y(s, t) = bsm(ks — cot), z(s,t) — 0, (57) 

for 0 < s < L where s is a parameter initialized as arclength. The wavenumber is k = 2tt/A for 



ihTvl 

-i_ 


3 4 5°1 2345 


(b) 


FIG. 6. Numerical results shown in the x-y plane for an undulating filament in a Brinkman fluid with 
7 = 10, b = 0.5, L = 5, A = 2, and uj = 27r. (a) The velocity field is shown along the structure, (b) The 
corresponding pressure map of the fluid domain. 


wavelength A, the bending amplitude is b, and uj is the constant angular speed. At any given time 
t > 0, the velocity of the flagellum is calculated by 


u x (s,t) — 0, u y (s, t) = — bu cos(fcs — cot), u z (s,t) = 0, (58) 

where u x ,u y and u z are the velocity components of x, y and z, respectively. Fig. [6](a) shows the 
sinusoidal swimmer with the velocity fields along the length of the swimmer in the x-y plane. The 
total velocity includes the velocity from the sinusoidal wave u s (x) = (u x . u y . u z ), the translation 
U 0 = (Uq, Uq , Uq), and the rotation of the filament f2 0 = (fig, fig, fig) as: 


v = u s (x) + u 0 + n 0 x x fc , 


(59) 


where V is defined similarly to Eq. (56) and for simplicity, we choose n — 1. Unless specified, the 
superscripts in translational and rotational velocity components are of the x, y and z components, 
not the partial derivatives. We note that f fc , U 0 and fl 0 are constants at each time point which can 
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be found by coupling Eq. (59) with the force-free and torque-free conditions. That is, 


V - U 0 - f2 0 x x fc = u s (x fe ), 


N 


5 > = °. 

k =1 
N 

^ ffc x x fc = 0 . 


( 60 ) 

(61) 

(62) 


k =1 


In Eq. (60), for each value of k, M £ is a 3A r x 3A r matrix while the coefficients for U 0 and fi 0 
will form 3 x 3N matrices. The coefficient matrices in Eq. and Eq. ( |62| ) are (3N + 6) x 3. 
To determine U 0 , f2 0 , and f, we solve Eq. 


|)-([62|). We can then compute pressure using the 
regularized version of Eq. ( |49| ). In Fig. [6jb), the pressure in the x-y plane is shown where we note 
larger variations in pressure close to the swimmer. 

The numerical results for the translational velocity U 0 will be used to compare to the asymp¬ 
totic swimming speed U ,^ derived in Eq. (|29j). The results are presented in Fig. [7} Hereinafter, 




(a) (b) 

FIG. 7. (a) The comparison between the asymptotic swimming speed (solid line) with the numerical data 
(marker points) for different amplitudes when the length of the swimmers are L = 5, L = 10, and L = 50 
and the permeability is 7 = 1. (b) The difference between the asymptotic values with the numerical results 
for different amplitudes, and different lengths for 7 = 1. The wavelength is taken to be A = 5 for all 
simulations. 


the wavelength is taken to be A = 5, t = 2 tt and u = 2n. We prescribe live different amplitudes 
for the simulations as b = 0.05, 0.075, 0.1, 0.125 and 0.15. We set the permeability 7 = 1 
and study the effects of swimmer length on swimming speed. We observe in Fig |7Ja) that the 
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numerical data (marker points) have good agreement with the asymptotic analysis (solid line) with 
a longer length (L = 50). We also plot the difference between the asymptotic values with the 
numerics for different lengths and different amplitudes when 7 = 1 as shown in Fig. [7j b). As 
the length increases, the differences decrease. This shows that finite-length swimmers will swim 
slower than the asymptotic predictions and this difference decreases for smaller amplitude b (with 
fixed oj and A). We also observe the same trend for different permeabilities 7 = 0.1 and 7 = 10 in 
Fig. [8]^a)-(b). We note that when looking closely at the error in Fig. |7Jb) for 7 = 1 and Fig. [H}a) 
for 7 = 0.1, the error is slightly larger for smaller permeability. Thus, the infinite-length cylinder 
swimming speed captures the swimming speed of a finite-length swimmer with more accuracy for 
larger permeability. 




(a) (b) 

FIG. 8 . Difference between analytical results and numerical data for different amplitudes and different 
lengths when (a) 7 = 0.1 and (b) 7 = 10 . 


In addition to the translational velocity, we can calculate the angular velocity for different 
parameters. We have five different amplitudes for the swimmer ranging from b = 0.05 to b = 0.15 
with five different permeabilities 7 = 0.01, 0.1, 1, 10 and 7 = 100. Fig. [9ja) shows that the 
angular velocity when L = 5 increases linearly as the amplitude increases. We capture the same 
behavior for longer swimmers (at L = 10 in Fig. |9}b) and L = 50 in Fig. [9}c)). We note that 
the angular velocity is much larger in the case of small length; in order for a swimmer of shorter 
length to achieve a prescribed amplitude, the angular velocity increases. 
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FIG. 9. Angular velocity for swimmers for five different permeabilities 7 = 0.01, 0.1, 1, 10, 100 with 
different amplitudes b = 0.05, 0.075, 0.1, 0.125, 0.15. (a) L = 5, (b) L = 10, (c) L = 50. 


2. Helical Bending Waves 


For this test case, we calculate the external torque exerted on the filament by the surrounding 
fluid. Consider the right-handed helix where the configuration is parameterized by the 3D space 
curve as 


N 



FIG. 10. Right-handed helix immersed in a Brinkman fluid with, n = 0.25, L = 20 and A = 5. The flow 
field is shown at z = 0 and z = 10. 


x(s, t) = ry cos(ks + ojt), y(s, t) — ry sin(fcs + out), z(s,t) = r 2 s + Ut, (63) 

for s, k, u defined as above, ry is the radius of the helix (or the amplitude), r 2 is a constant defined 
as r 2 = cos 9 = y/l — k 2 r\ where 9 is the pitch angle, and U is the constant propulsion velocity. 
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The prescribed helical configuration gives the velocity of the helix as 

u(s, t) = — r 1 cusin(/cs + cut), v(s,t) — ricu cos(ks + cut), w(s,t)=U. (64) 

The torque is calculated as [[480 

T = J SkX ds, (65) 

where T is the helix (centerline of the flagella) and gy is the surface force (traction) applied on the 
filament. The torque is then numerically approximated by 


N 


T = ^2 (g k X Xfc) As, 


( 66 ) 


k =1 


which we will compare with the analytical solution T,' x in (42). 

In Fig. |TT] we present results for the helical bending case for amplitudes (ry = b) in the range 
of 0.05 to 0.15 and set the wavelength at A = 5. For a filament of radius a = 0.05 and permeability 
7 = 0.01, we observe an increase in torque as amplitude increases in both the asymptotic anal¬ 
ysis and the computational solutions. We note that similar to the swimming speeds for the case 
of planar bending, the asymptotics greatly overestimate the torque for shorter length filaments. 
Similarly, for the cases of 7 = 1 and 7 = 10 (shown in Fig. |TT} b)-(c)), the analytical results for an 
infinite-length spiral cylinder overestimate the torque for the finite-length spiral filament. The dif¬ 
ference between the numerical data and analytical solutions is also plotted in Fig.|TT|d) for 7 = 10 
to show that for each fixed amplitude 77, as the length increases, the difference decreases. We note 
that previous computational studies using the MRB have observed that the optimal numerical reg¬ 
ularization parameter e varies for each 7 and can be sensitive for torque calculations li44l . We also 
observed this sensitivity and decreased the regularization parameter as permeability increased to 
report the best fit with the asymptotics at longer lengths (values given in Figure caption). 


VIII. DISCUSSION 


In this paper, we have analyzed an infinite-length cylinder undergoing periodic bending in a 
fluid governed by the Brinkman equation, which is a model for flow through a porous medium. 
Motivated by organisms that exhibit undulatory locomotion such as spermatozoa, we focus on the 
case where the radius of the cylinder is small in comparison to the fiber spacing. We find that 
propulsion in the case of planar and spiral bending is enhanced with larger fluid resistance for spe¬ 


cific combinations of wavenumber, cylinder thickness, and permeability. In Section VI we show 
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FIG. 11. The comparison between the asymptotic velocity (dash-line) with the numerical data (marker 
points) for different amplitudes and different length, L = 10, 20, 30, 40 and L = 50 with a = 0.05 when 
(a) 7 = 0.01 and s = 0.01, (b) 7 = 1 and e = 0.0055, and (c) 7 = 10 and e = 0.005. (d) The difference 
between the asymptotic values with the numerical results for different amplitudes, and different lengths for 
7 = 10. 


that for a sufficiently small volume fraction of fibers, mammalian spermatozoa will observe an en¬ 
hancement in swimming speed. Our calculations show that the mesh spacing is several orders of 
magnitude larger than the cylinder thickness, allowing room for the swimmer to navigate between 
stationary fibers. This analysis provides insight into the sperm thickness and wavelength that we 
observe in nature; perhaps they have been optimized to provide enhanced swimming speeds in 
oviductal fluids. 

The observed enhancement in swimming speed for the infinite-length cylinder with planar 
bending is similar to the case for a swimming sheet in a Brinkman fluid If35l1 . We note that in both 
the 2D and 3D cases, as the resistance is reduced to zero, the corresponding swimming speeds in 


25 





















a Stokes fluid are recovered. For a fixed amplitude, wavenumber, and cylinder thickness, the ratio 
between the swimming speeds of the infinite-length cylinder and sheet is approximately 0.8 in a 


Stokes fluid (using Eq. (J33])). We observe that the ratio of the asymptotic swimming speeds of the 
3D infinite-length cylinder and 2D sheet in a Brinkman fluid vary greatly as the scaled resistance 
increases, decreasing from 0.8 to 0.1. Thus, as the scaled resistance increases, potential rotational 
effects may play a large role in decreasing the swimming speed of the 3D infinite-length cylinder 
(in comparison to the 2D sheet). This highlights the importance of the analysis presented here for 
the 3D infinite-length swimmer, especially when trying to understand the role of larger resistance 
on swimming speeds. 

Relative to the Stokes case, the infinite-length sheet and cylinder swim slower in a viscoelastic 
fluid Il24l - t26ll . However, the results reported here for an infinite-length cylinder and previous 
work for the sheet 051 show that added fluid resistance enhances swimming speeds relative to 
the Stokes case. Thus, a potential enhancement in swimming speed can be observed when a low 
volume fraction of obstructions do not have a frequency dependent response or when the polymer 
relaxation time is fast. In contrast to a viscoelastic fluid, the Brinkman fluid model assumes that 
the fibers or polymers in the fluid are stationary. 

Through a detailed mathematical analysis, we have derived the work for planar waves. We 
observed that larger resistance (smaller permeability) results in a large increase in work. This 


increase in work will occur when the cylinder thickness satisfies the inequality in Eq. ( |45] ). We 
note that in the asymptotic derivation, we have assumed prescribed kinematics. Thus, as the 
resistance increases, it requires more work to maintain planar bending with the same amplitude and 
wavenumber. When building artificial microswimmers in a porous medium, one must consider the 
amount of energy required to have the swimmer bend f[49ll . This could be a constraint on reaching 
higher swimming speeds in fluids with larger resistance. 

We have compared our asymptotic solutions to computations of finite-length swimmers with 
prescribed kinematics using the method of regularized Brinkmanlets. For cylinders of sufficient 
length, the asymptotic swimming speeds match well with the computations. We note that the 
asymptotic analysis is able to capture the trends of swimming speed in terms of the dependence 
on permeability and amplitude. However, it overestimates the swimming speed for shorter length 
filaments. This is important to consider when using asymptotic swimming speeds to make predic¬ 
tions of the behavior of finite-length swimmers. Additionally, we have observed that the analytical 
results overestimate the torque for a finite-length filament with helical propagating wave. 
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For prescribed kinematics, we note that the asymptotic and computational swimming speeds 
calculated increase as amplitude and resistance increase. The asymptotic analysis presented here 
provides the swimming speed given that a swimmer could attain the given prescribed kinematics in 
a fluid with permeability 7 . Previous studies have observed non-monotonic changes in swimming 
speed for finite-length swimmers with increasing fluid resistance for planar swimmers, where the 
achieved amplitude of bending is an emergent property of the fluid-structure interaction 144, ;501 . 
In these studies, finite-length swimmers were not able to achieve large amplitude bending as the 
permeability is decreased. In addition, experimental studies have shown that the emergent wave¬ 
forms and swimming speeds will depend strongly on the fluid environment [T3J EH]. Thus, it 
is important to put the asymptotic results in the context of finite-length swimmers where certain 
ranges of bending kinematics are not observed in gels or fluids with small volume fractions of 
fibers. 

In this computational study of a finite-length filament undergoing periodic lateral bending, we 
observed a large increase in angular velocity as the swimmer length decreases. Additionally, angu¬ 
lar velocity increased linearly as amplitude increased for a fixed beat frequency. Sperm cells have 
been observed to ‘roll’ as they swim (simultaneous rotation of the sperm cell body and flagellum) 
[13. ED. Specifically, human sperm were found to increase rolling from 1.5 Hz to 10 Hz and 
decrease amplitude as the viscosity of methylcellulose solutions was decreased lfT3ll . In our com¬ 
putational study, angular rotation (rolling) varies linearly with amplitude and is much smaller than 
the experimental data. However, we are not accounting for the dynamics of a cell body and have 
prescribed kinematics. It will be interesting to study three-dimensional computational models of 
finite-length swimmers with cell bodies and emergent kinematics in the future to fully understand 
swimming speed and angular velocity as a function of the permeability. 


IX. APPENDIX 

A. Derivation of surface cylinder location for planar bending waves 

The surface points of the cylinder are defined by Eq. ([3]). Using Fig.|T| we observe that O' = 0+£ 
and use this to rewrite the equation as: 

r 2 = a 2 + b 2 sin 2 s + 2 ab sin s cos (0 + £) 

= a 2 + b 2 sin 2 s + 2ab sin ^(cos# cos£ — sin 6 sin£) 
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, . ^ b sins sin 9 

where sm £ = -. For small we have: 


r 2 = a 2 + hi 2 sin 2 s + 2 ab sin s cos 9 — 2 b 2 sin 2 s sin 2 9 
= a 2 + 2 ab sin s cos 9 + b 2 sin 2 s(l — 2 sin 2 6 ) 

= a 2 + 2 ab sin s cos 6 + b 2 sin 2 s cos 2 9 — b 2 sin 2 s sin 2 6 
= (a + b sin s cos 0) 2 — b 2 sin 2 s sin 2 9 
= (a + b sin s cos 9) 2 


r = (a + b sin s cos 9) y 1 — 

We can then arrive at the final equation for r: 


1 - 

b 2 sin 2 s sin 2 9 

(.a + b sins cos 9) 2 


/ 

b 2 sin 2 s sin 2 9 


( 1 a + b sin s cos 9 ) 2 


r = a 

1 H— sin s cos 9 + 0 

m 2 ’ 


a 

W 


(67) 


B. Asymptotics 


We wish to calculate the swimming speed of the infinite cylinder. The velocity components are 
expanded up to the second order about C = Ci = ka: 

u = (w)^=f 1 + bk cos 9 sins (V)^^ + • • • , 

v = (v)c=Ci + bk cos 9 sin + • • • , (68) 

w = (w)^=^ + bk cos 9 sin s{w')^= Ci > 

where Eq. Q is used to rewrite (( — Additionally, the velocity components u, v, and w are 
expanded in the powers of 6/a, 


u = u + • • • , v — 7 / 1 ) + r/ 2) + • • • , w — w ^ + w^ + • • • . (69) 

Substituting Eq. ( [69] ) into Eqs. ( [ 68 ] ) and ( fTO] ): 

u — u 1 cos 9 cos s + u ^ | + bk cos 9 sin s cos 9 cos s(m / )^ =( ^ 1 , (70) 

v = V\ sin 9 cos s + v ^ | + bk cos 9 sin s sin 9 cos s(v')^ = Cd (71) 

w — wi cos 9 sin s + + bk cos 9 sin s cos 9 sin s(w')^= Ci- (72) 


By matching the above expansions with the boundary conditions in Eq. Q, we can determine the 
constant coefficients A mn , B rnn and C mn for each order of the expansion. 
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Specifically, the second order expansion is: 


uW 

w ^ 


—bk cos 9 sin s cos 9 cos 

—bk cos 9 sin s sin 9 cos s{v')^=^ 

—bk cos 6 sin s cos 9 sin s(w')q = ^ 


— -6/c('ti / 1 )^ = ^ 1 (sin2s + 008 26*811125), (73) 

—jbk(v\)^=^ sin 29 sin 2s, (74) 

— -bk(w' 1 )^ 1 (l — cos 2 s + cos 2 9 — cos 29 cos 2s). 

(75) 


The coefficients of the velocity in the second order expansion can be evaluated as: 

“0 2 = ~^k(u[) c=Ci , u 22 = -|^K) c=Cl , 

V22 = ~^bk(v[) c= Cl , 

1 11 

w 0 2 = -&fcK) c=Cl , w 20 = --bk{w[) c=(:i , w 22 = -bk(w[) c=<1 . 

The cylinder will move at a speed of (Too with respect to the fluid at infinity which is also the term 
that balances the constant expression in w^ 2 \ Then, 

(Too ~ l)k (w ] ) =fca • 

C. Derivation of surface cylinder location for spiral bending waves 

As shown in Fig. |4j the time-dependent contour caused by the spiral bending wave is defined 
as 


or equivalently, 


a 2 = r 2 + b 2 — 2br cos(6* — s) 

= [r — 6cos(6* — s)] 2 + b 2 sin 2 (9 — s ), 


r = bcos(9 — s) + \j a 2 — b 2 sin 2 (9 — s ), 


r = a 


1 H— cos(6* — s) + O ( - 


29 





D. Analysis of the Asymptotic Functions 


We look more closely at the behavior of the velocity in Eq. (29) and the work done in Eq. (35). 
Rewriting in terms of the scaled resistance ft = a/k. 


m = 


u n 


x ° (ci)- K^ +1 ) iog(i+/?2) 

< 17 ^ = Ao(G) + + 


gift) = 


W 


^ b2k2u2 mi) + U±-i)io g a + P 2 ) 


(76) 


(77) 


The two functions are plotted in Fig. 12 Using the condition in (44), f(ft) and g(ft) are positive 
functions and f(ft) is bounded by 1. The first derivatives of f(ft) and g(ft) with respect to ft are 


m = 

9 {8) = 


2ft 2 K 0 (Cr) 


— 1 + log(l + ft 2 ) + -g 2 l°g(l + ft 2 ) 


-o(l + ^) l°g 2 (l + /3 2 ) 


(1 + ft 2 


1 / 1 


2 \ft 2 

(1 + ft 2 ) log(l + ft 2 ) + /3 2 (—1 + ft 4 ) 


I<o(Ci) + - -,-1 Uog(l + 8 2 ) 


(1 + ft 2 ) 


MCi) + \(jp ~ 1 ) log (! + ft 2 ) 


1 2 ‘ 


, (78) 


(79) 


We observe that all terms in the denominator and the numerator of g'{ft) are always positive for 
all ft which implies gi ft) is an increasing function. On the other hand, the function inside the curly 
bracket of f'(ft) is positive when 

(1 + 1/ ft 2 ) log 2 (l + ft 2 ) 


K 0 (Ci) > 


2[—1 + log(l + ft 2 ) + ( 1/ft 2 ) log(l + ft 2 )]' 


( 80 ) 


In other words, f(ft) is an increasing function when it satisfies the condition in (80). We note that 
the expression —1 + log(l + ft 2 ) + (1/ft 2 ) log(l + ft 2 ) is always positive. The Taylor expansions 
of f(ft) and g(ft) about ft -C 1 are as follows: 

1 




gift) 


Ko(Ci) + 1/2 


+ 0(ft 2 ). 


( 81 ) 


This shows that when ft is small, f(ft) > g(ft) as in Fig. 12 When ft is large, we can expand the 
two functions in terms of the Puiseux series as: 


m 

gift) 


i + 


2 log(l/ ft) 


P [K 0 (Ci) + log(l/ft)] 




Kq((i) + log(l/ ft) 


+ 0i p 


(82) 

(83) 
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FIG. 12. The plot of the functions / in Eq. ( |76| ) and g in Eq. ( f77] >. 

Clearly, f(/3) is bounded by 1 when (3 is large while g(/3) is unbounded. The two formulations 
above give insight as to why a decrease in permeability 7 causes a small increase in swimming 
speed and a large increase on the rate of work done. 
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